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Abstract 

We discuss the stabilization of the conformal factor by higher derivative terms 
in a conformally reduced R + Euclidean gravity theory. The flat spacetime is 
unstable towards the condensation of modes with nonzero momentum, and they 
"condense" in a modulated phase above a critical value of the coupling /3 of the 
R^ term. By employing a combination of variational, numerical and lattice meth- 
ods we show that in the semiclassical limit the corresponding functional integral 
is dominated by a single nonlinear plane wave of frequency ~ l/y^^pi. We argue 
that the ground state of the theory is characterized by a spontaneous breaking of 
translational invariance at Planckian scales. 



1 Introduction 



It is not uncommon in euclidean field theories that spatially inhomogeneous, i. e. non- 
constant field configurations give rise to a lower value of the action functional than homo- 
geneous ones. As a result, at least at a semiclassical level, inhomogeneous configurations 
are likely to dominate the functional integral and to determine the quantum vacuum 
state, |0). In this situation the main properties of the true quantum state often can be 
understood by a saddle point expansion of the quantum fiuctuations {i.e. , the integration 
variable of the functional integral) about a specific set of field configurations; the lat- 
ter have a position-dependent, translation invariance-breaking value of the field variable, 
and they are degenerate with respect to their value of the action. In the quantum theory 
which is approximated in this way, this "condensation" of spatially inhomogeneous modes 
contributes to certain expectation values (0|(9|0) 7^ which are sensitive to the nonvan- 
ishing kinetic energy of the configurations dominating the functional integral. Here O is 
a scalar operator constructed from the fundamental fields; for instance O = d^(f)d^(j) in a 
scalar model, 01 O = ti^F^^F^'^) in Yang-Mills theory. Such contributions are sometimes 
referred to as "kinetic condensates" [Ij. They are to be distinguished from the more famil- 
iar translational invariant "potential condensates" which underlie the conventional Higgs 
mechanism; there the functional integral is dominated by nonzero but spatially constant 
scalar field configurations. 

Let us denote the fundamental fields collectively by $, the classical action by 
its stationary point(s) by $0; ^•C-|f['^'o] = 0, and its global minimum by $mm- Then 
the inverse propagator which governs small fiuctuations about a configuration $0 is given 
by the Hessian S^^^f^] = 5S'[$]/(5$5$ evaluated at $ = $0. If the operator ^(^^[^o] 
has negative eigenvalues then there exist directions in the field space along which the 
action can be lowered. This leads to a run away behavior of certain fiuctuation modes 
whose excitation (condensation) brings the field close to the global minimum of the action, 
$min. Fluctuations about ^min instead should all be stable, i.e. the Hessian at the global 
minimum, S''-^^ ['^'min] , has only non- negative eigenvalues. 

When one uses stability criteria in order to judge whether a theory is physically 
acceptable, one may not confuse the Hessian at some arbitrary solution of the field equa- 
tion, 5'*^^^[$o]) with the one at the absolute minimum, 5**^^^ [$min] • Only the latter must 
have a positive semidefinite spectrum. In the example where $ is the Higgs field, the 
homogeneous configurations $o(a;) = and $min = say, are both stationary points, 
corresponding to the "false" and the "true" vacuum respectively. But clearly only the 
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S'*^^'' [$min] is positive; the Hessian S'*^^^[$o] has a negative eigenvalue corresponding pre- 
cisely to homogeneous fluctuations which tend to drive the fleld from $ = to $ = f . 

In theories with instabilities of the kinetic type the situation is conceptually similar, 
albeit more complicated technically. Here the transition from a false to the true vacuum 
involves not just a shift of the fleld by the vacuum expectation value f , as in the Higgs 
case, but rather an expansion about a position dependent fleld conflguration. 

A well-known example of a kinetic condensate is the gluon condensate in Quantum 
Chromo-Dynamics (QCD). There, the classical action functional |J df'xii {F^yF^") is 
minimized by gauge fleld conflgurations with F^^^ = 0, but the effective action has its 
minimum at F^i, ^ 0. An early attempt at flnding its global minimum is the Savvidy 
vacuum [2], the approximation of a covariantly constant color magnetic fleld. While its 
action is indeed lower than that of the naive vacuum with F^,^ = 0, it turned out unstable 
in the infrared (IR), and it has been argued that the true vacuum should be spatially 
inhomogeneous. The complexity of the QCD vacuum state is reflected by nonperturbative 
contributions to (0|tr (F^j.F^'^) |0) and similar expectation values of more complicated 
gauge and Lorentz-invariant operators. 

Typical examples from statistical physics are materials described by a Landau free 
energy functional with a Lifshitz point |3] and, among the classes of magnetic materials, 
the antiferromagnetic ones. These latter display a microscopic order characterized by anti- 
parallel sublattices of spins, so that their global moments are exactly equal but opposite. 
Below the Neel temperature the susceptibility obeys the Curie- Weiss law for paramagnets 
but with a negative exchange interaction. The antiferromagnetic order can then show 
up as spatial inhomogeneity [1]. Similarly, in a superconductor, the superconducting 
order parameter displays a modulated phase in the presence of a strong magnetic fleld, 
the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase [5]. (The antiferromagnetic phase in 
scalar fleld theories has also been described in [6].) 

The model we are interested in here is euclidean quantum gravity [7\. Besides 
perturbative nonrenormalizability, the other property the Einstein- Hilbert action 



is notorious for is its unboundedness below. Indeed, if we make the conformal factor of 
the metric explicit by setting g^i, = exp(2cr)^^,^ we obtain for the action 
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and its value can become arbitrarily negative when cr{x) varies rapidly. 
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Figure 1. The classical inverse propagator 0,{p'^) in eq. (jl.4p and the effective inverse 
propagator Oefr(p^) according to the toy model of eq. ljl.Sp . (Taken from [1].) 

However, it is easy to modify the action in such a way that it becomes impossible 
to lower the action below all bounds by exciting the conformal factor; hereby all the 
successfully tested predictions of classical general relativity are retained. The simplest 
way consists in adding a term proportional to the square of the Ricci scalar, with a 
positive coefficient (3 > : 

S[g] = SEH[g] + (3 j d^x^R" (1.3) 

The term involves 4 derivatives of metric. In an expansion about flat space, the mod- 
ified action implies an inverse propagator for the conformal factor which, schematically, 
is of the form fi(— □) = □ + □□/(2M^) or, in momentum space, 

^ip') = + 1^ (1-4) 

Here the —p^ and terms are due to the Einstein-Hilbert and term, respectively, 
and M oc mpi/ is a constant of the order of the Planck mass. For momenta much 
smaller than M the —p^ term dominates, and the field can lower its action by exciting 
such modes. If ^ M^, on the other hand, the positive -|-p^/2M^ term dominates and it 
"costs" action to excite the corresponding modes. For euclidean momenta = p^p^ > 
the kinetic operator fi(p^) is positive for p^ > 2M^; it has a minimum at p"^ = where 
it assumes the value f2(M^) = —^M^. So, according to the modified action, the kinetic 
energy of the conformal factor is bounded below. 

The important point is that the field configuration corresponding to fiat space {g^^ = 
5^v) cannot be at the absolute minimum of S. In fact, fiat space plays the role of a false 
vacuum here. While in absence of a cosmological constant it is a stationary point of S, 
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the Hessian S^'^^[g^„ = 5^^^] has negative eigenvalues which correspond precisely to the 
fluctuation modes of the conformal factor with fi(p^) < 0. This is the typical symptom 
of a "kinetic condensate" that wants to form in order to lower the field's value of the 
action functional. So in order to analyze both the classical and quantum properties of the 
theory based upon the modified action fll.3p it is important to have some understanding 
its minimum action configuration(s). 

Before we turn to this problem several remarks are in order. 

(A) The perturbative quantization of i?^-gravity in an expansion about fiat space results 
in a power counting renormalizable, though non- unitary theory [8] . This well-known fact 
by no means rules out the existence of a bona fide quantum field theory based upon the 
action. It simply says that flat space is not the true ground state but only a false vacuum 
of the theory; the negative-norm states ("ghosts") one encounters expanding about this 
false vacuum are a reflection of unstable eigenmode of 5*^^^ [g^^ = 5fj_u] which tend to grow 
when the system is heading for its true ground state. The expansion about the true 
vacuum can very well be stable and unitary. 

(B) The R^- addition to the Einstein-Hilbert term is far from unique. In this paper 
we study it as the simplest example of a theory in which the kinetic energy due to the 
conformal factor is bounded below. The basic condensation mechanism should be similar 
in all models where ^^(j9^) = — + /(p^), with f{p^) > a monotonically increasing 
function for — )■ oo. 

(C) Besides the conformal factor problem, the other key issue in trying to construct a 
quantum field theory of the metric is the ultraviolet (UV) renormalization. Even though 
this is not directly relevant for the investigations in this paper, we assume that the ultimate 
(stable) theory can be constructed along the lines of the Asymptotic Safety program [9] 
based upon the gravitational average action [10] . The basic idea is to take the limit of an 
infinite UV cutoff at a non-Gaussian renormalization group fixed point. The resulting field 
theory can then be predictive and well behaved on all scales. By now there is significant 
evidence for the existence of a suitable fixed point [TTH39] . More than that, it has also 
been found [14] that the corresponding bare action seems to contain an i?^-term with a 
positive coefficient precisely what is needed for the envisaged condensation mechanism. 

(D) In ref. [I] the condensation of spatially inhomogeneous modes has been studied in 
detail within a scalar toy model which mimics certain features of R + R^ gravity. It 
consists of a complex scalar field X; on c?- dimensional fiat euclidean space, governed by 
the action 




(1.5) 
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The kinetic operator is the same as in Eq.f ll.4p . so exciting the x-iiiodes with momenta 
in the interval G [0, 2M^] lowers the action. It was shown that the global minimum of 
the functional fll.Sp is given by the family of plane waves 



labeled by a unit vector n G S'^^^ and a phase angle a. The classical vacuum manifold is 
S'^^^ X 5^ therefore. 

In ref. flj the dressed inverse propagator Vtcs_{p'^) which appears in the model's 
effective action functional T[x\ = J x*^efr(— n)x + ■ ■ • was computed by a saddle point 
expansion about the configurations (11. 6p which involved an integration over and a. 
The result for fieff(p^) is depicted in FigJU It shows quite nicely the dynamical self- 
stabilization of this theory by the condensation of spatially inhomogeneous modes: due 
to the renormalization effects, the kinetic term has become positive semidefinite; for all 
modes with ^ it "costs" energy (action, actually) to excite them. Only the modes 
with = can be excited "for free" ; this indicates that those modes might be unstable 
towards condensation. 

The analysis of this toy model could be carried through in a rather complete way, 
including the calculation of its effective average action Tk [x] interpolating between S and 
r. This was possible because of certain algebraic simplification arising from the very 
special form of the classical action, and in particular since x was taken to be a complex 
field. 

The present paper is intended to be a first step towards a similar analysis for a 
realistic gravity action, with a real conformal factor in particular. 

In the following we shall search for the global minimum of the R + B? action func- 
tional (with a cosmological constant included) and try to establish its modulated nature 
for appropriate values of the parameters in the action. This is a very hard problem 
which cannot be solved exactly. In order to make progress, we restrict the domain of the 
functional to conformally flat metrics, z.e. metrics conformal to M^. We shall then use 
variational and numerical methods to find approximations to the global minimum in this 
subspace. 

We believe that, as far as a possible dynamical resolution of the conformal factor 
problem is concerned, this restriction still contains the essential physics. In fact, the other 
( "transverse" ) metric degrees of freedom which we discard show no comparable instability 
and have no obvious reason to condense. 

We shall be particularly interested in finding periodic conformal factors which par- 
tition the spacetime into an array of elementary cells. If they exist, they could possibly 



X(x;n,a) 
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exp {iMn^x^ + ia) 



(1.6) 
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serve as a classical approximation of quantum Minkowski space, or its euclidean coun- 
terpart, in the following sense: On microscopic, typically Planckian scales the metric is 
violently oscillating, but upon averaging it over a periodicity volume it becomes perfectly 
fiat, ie. after a purely classical "coarse graining" one has {g^y{x)) = 5i^,y. 

The remaining sections of this paper are organized as follows. In Section 2 we dis- 
cuss the general properties of the conformally reduced R + action. In Section 3 we 
find its explicit global minimum for a special point in parameter space which saturates 
a kind of Bogomolny bound. In Section 4 we employ variational techniques to obtain an 
approximation to the global minimum when the cosmological constant is zero, and in Sec- 
tion 5 we use numerical techniques to illustrate the impact which a nonzero cosmological 
constant has on them. Section 6 contains the conclusions. 

2 The conformal sector of R + gravity 

In the following we consider the euclidean action functional 

S{g,u] = j rf'xy^|^(-i? + 2A) + /3i?2| (2.1) 

the cosmological constant A and the dimensionless parameter /3 are assumed positive 
throughout. We shall see that /3 > guarantees that S is bounded below. All stationary 
points of (12. ip . its global minimum in particular, satisfy 

[G^, + Ag^,] + P [-{G^, + R^,)R + 2D^D,R - 2g^,D^R] = (2.2) 

with the Einstein tensor G^y = R^y — \g^uR- Contracting (12. 2p with g^^" yields 

- 4A + (967rG/3) D'^R = (2.3) 

A special feature of the R^ action (12. ip is that every stationary point of the Einstein- 
Hilbert action (/3 = 0) continues to be a stationary point when the [iR^-term is added. In 
fact, it is easy to check that metrics satisfying the ordinary Einstein equations 

G^y = -Kg^y ^ R^y = Ag^y, R = AA (2.4) 

automatically solve eq. (l2.2p also for (3^0. 

Mostly we shall be interested in those stationary points of S which are not already 
stationary points in absence of the stabilizing R^ term. As a consequence, (12. 4p will not 
be satisfied in these cases, and eq. (12.31) tells us that the curvature scalar is non-constant 
then. More precisely, R{x) may not be a harmonic function, D^R ^ 0. 
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By completing the square in eq. fl2.ip we can rewrite the action in the form 

7 

(327rG')2/3 



S[g] = I 



' + TT^TT-^ } (2-5) 



with the abbreviation 

7 = 1287r /3GA - 1. (2.6) 

The representation fl2.5p is vahd if I3GA 7^ 0, or 7 7^ —1. Since, by assumption, /3 > it 
imphes a lower bound on the value of the actionl^ 

S[g] > f d'x^ (2.7) 



(327rG)2/3 

Obviously the functional S is positive definite if 7 > 0, or 

1287rG'A > 1. (2.8) 

If 7 = 0, that is, for the special parameters values for which 128tt/3GA = 1, the 
minimum value of the action is at 5 = 0, and the bound is saturated by metrics satisfying 
R = 4A. While this is precisely the contracted form of the ordinary Einstein equations, 
here R = AA has the interpretation of a kind of "Bogomolny equation" . 

In the case 7 < the functional S can assume negative values, but it is still bounded 
below. 

Rather than trying to find the global minimum of S[gf^,y] defined over the space of 
all metrics we shall be more modest here and only analyze the action restricted to the 
conformally flat metrics 

= (2-9) 

Writing 5'[0] = S[gfj_,, = (p'^S^i^] for the functional depending on the conformal factor we 
obtain from (12.11) : 



Sl^]= li'.{^^n^^—^' + ,6p[^^) } (2.10) 

where we used that y/g = 0^ and R = — 60^^n0 with □ = for metrics of the 

form g^y = (t)^5^u- If /? = the restricted functional (I2.10p has the appearance of a scalar 
</)^-action with a "wrong sign" kinetic term. 

The main topic of he present paper is the investigation of the global minimum action 
configuration(s) 0min of the restricted functional S'[0], Eq. (l2.10p . It is plausible to assume 



"'^We consider compact manifolds without boundary or, in the infinite but periodic case, we refer the 
action to a single periodicity volume. Hence J d*Xy^ is always finite. 
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Figure 2. Schematic view of the restricted functional space studied in the paper. The 
paraboloid represents the action for the complete set of degrees of freedom, its global mini- 
mum at the origin thus satisfies the field equation for the R + B? model. The plane of frozen 
vector and tensor degrees of freedom cuts out instead a parabola (in bold in the figure) 
which represents the restricted functional space of the conformally reduced theory. (The 
axis with the arrow is the direction of the conformal field configurations.) 



that its essential qualitative features, in particular the existence of a modulated phase 
for certain parameter values, will be shared by the true minimum i.e. that of S[g^y\ 
defined for all, not necessarily conformally fiat metrics. 

While we hope that g"^^^ is to some extent similar to the conformally flat metric 
with the lowest action, 

<^>lM8,u (2.11) 



conf— min 



we emphasize that g^^ 
(12. lip is found from 



conf —min 



is not a stationary point of S[g^y\ in general. The metric 
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6(j){x 



-S[4>] 



6(j){x 



-S[g, 



fJ,U 



(2.12) 



</'=0min 



while c/™"" satisfies (^^, i.e. 
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■S[g] 







(2.13) 



Sg^^uix) 

Fig. ([2]) is a graphical illustration of the fact that (12. lip with fl2.12p does not imply ( 12.131) 

^con 



in general. In fact, the metric g'^J^^ sits to the minimum of the parabola cut out by 
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the plane, but it is not a stationary point of the unrestricted functional represented by 
the paraboloid. 

(A) Before turning to the actual minimization problem a remark concerning dimensions 
might be in order. Throughout this paper, the coordinates denoted x'^ are assumed 
dimensionless. Hence, in this system of coordinates, all metric components g^^ have mass 
dimension —2, for ds^ = g^i,dx^dx^ has [ds^] = —2 always. As a result, the conformal 
factor introduced as in (12. 9 p has the dimension of a length, [0(a;)] = —1. 

(B) Sometimes it is advantageous to introduce the conformal factor in a slightly different 
way so that the second order kinetic term in S has a standard normalization (up to its 
sign). Writing 



g^,{x) = ^(AnG) x'ix) V, ^-e. 0' = ^(47rG) (2.14) 



the action for x reads 



Six] = / rf^x^ - XOX + 36/3 j +-x'\ (2.15) 



36/3 / d'xxH^ + T^y + I d'xx' (2.16) 



V x^ 14:4/3 J 576^ 



with the parameters 



u = —GA, 7 = 1287r PGA - 1 = 24 m/3 - 1 (2.17) 
3 

As [G] = —2, [A] = +2 the parameter u is dimensionless, the same is true for the new field 
variables: [x] = 0. Since ip\ = \fG = rripl is the Planck length, we see that x measures 
proper distances ds"^ = ^{x{x)ipiy S^^^dx'^dx'^ in units of £pi. 

(C) The above convention of dimensionless coordinates is most convenient when dealing 
with curvilinear coordinates, in curved space in particular. It is, however, not the con- 
vention usually adopted in ordinary quantum field theory on Minkowski space. There 
one prefers using cartesian coordinates with respect to which the metric components are 
normalized to ±1. Since ds"^ is still required to have dimension —2, those coordinates 
necessarily have the dimension of a length. 

The transition from our above conventions to this kind of dimensionful coordinates, 
henceforth denoted x'^, is achieved by the rescaling 

[x^^] = -l (2.18) 
In terms of the x's the line element implied by a function x(3^") reads 




ds^ = -,x\x- = ^^^.) S.^dx^dx'^ (2.19) 
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Here c > is an arbitrary dimensionless constant which we may use to change the absolute 
normahzation of the x's. For instance, if happens to be a periodic function we might 
choose (? = (x^(')) to be the average value of taken over one periodicity volume 
(elementary cell). As a result, the line element if averaged correspondingly is precisely 
that of flat euclidean space expressed in terms of standard cartesian coordinates: 

{ds^) = S^^dx^'dx" (2.20) 

Later on we shall indeed find minimum action configurations x{^) which have a Planck 
scale periodicity and average to fiat space on macroscopic scales. 

(D) Varying the action S'[x], eq.f l2.15p . with respect to x yields the following 4th-order 
partial differential equation for its stationary points: 



o 



(2.21) 



Besides the global minimum of S\x\ which we are after, this equation is satisfied also by 
many local minima and maxima and stationary points of mixed type. In particular all 
solutions xo of the equation (12.211) with /3 = are solutions of the full equation with 
/3 7^ 0, too. In fact, we see immediately that 

□Xo+gXo = (2.22) 

implies /3[n(xo ^Dxo) ^ Xo ^(OXo)^] = 0. Thus the solutions of this somewhat trivial type 
are insensitive to the value of /3. 

(E) In general the partial differential equation f l2.2ip . equipped with appropriate bound- 
ary conditions, is hard to solve. In Section 5 we shall analyze it using numerical methods. 
For that purpose, still another parametrization of the conformal factor, namely in terms 
of an exponential turned out advantageous: 

(7^, = ^(47rG')e2-(-)<5^„ i.e. ^(x) = e'^(^). (2.23) 

The corresponding action reads 

S\g\ = j dS {^e2'^(na + d^ad^a) + 36/3 (Da + d^ad^af + ^e^'^} (2.24) 

and the condition for stationarity assumes the form 

e'^fDa + ^^a^^'a\ + g e^" + 72/3 [nDa + 2{d^d,a){d^d'' a) - 2{Ua){Ua) 

-2{d^(y){d^'a)U(y - 4(9^a)(9,a)(9^9V)] = 0. (2.25) 
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In the rest of this paper we shall approach the minimization problem of S by three 
different methods: In Section 3 we study the special case 7 = by means of the "Bogo- 
molny equation", in Section 4 we use a variational method, and in Section 5 numerical 
techniques. 



3 Special case 7 = 0: The Bogomolny equation 

In this section we assume that the parameter combination 7 = 128tt/3GA — 1 assumes the 
special value 7 = 0. Then Eq. (12.51) and Eq. fl2.7p imply that > and equality holds 

for metrics with 

R{g) = 4A (3.1) 

We shall refer to (13. ip as the "Bogomolny equation". It happens to coincide with the 
contraction of the ordinary {i.e. second order) Einstein equation. If 7 = 0, all solutions of 
(13. ip saturate the lower bound S[gf^u] = 0. 

Here we shall consider only conformally flat solutions to Eq. (l3.ip . Writing the metric 
as in Eq. (l2.14p and inserting it into (13. ip we are thus led to solve the Yamabe problem 

nx + -x' = (3.2) 

This equation has the same mathematical structure as Eq. (12.220 . its interpretation is 
somewhat different though. To get (12.221) from (12.211) we had set /3 = 0, but (13.21) refers 
to the specific nonzero (3 for which 128n f3GA = 1. (Recall also that u = ^GA whence 
u/6 = 1/(144/3) if 7 = 0. Therefore Eq. (l3.2p coincides with the Bogomolny equation one 
reads off from the reduced action (I2.15p .) 

The only known periodic solutions to equation (13. 2p are the euclidean analogs of 
traveling plane waves. If x depends only on one of the coordinates, = x, say, Eq. (l3.2p 
can be interpreted as the Newtonian equation of motion of a particle moving along the 
"X-axis" under the influence of a potential V{x) = • 

x"{x) = -^Vix) (3.3) 
dx 

Here the prime denotes a derivative with respect to x which plays the role of time. Since 
the cosmological constant is assumed positive, u and V are positive too, so that the 
solutions to (13. 3p correspond to anharmonic oscillations. They can be found explicitly in 
terms of Jacobi elliptic functions [H]. 

The general plane wave type solution of (13. 2 p involves two free constants, x and a. 
It reads 



^(^x) = xsn{x\l ^rif.x^' + a;i) (3.4) 
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Here sn denotes the sinus amplitudinis with purely imaginary modulus k = i, and n'^ is 
an arbitrary unit vector, 5^yn^n'^ = 1. Obviously the solutions (13 ■4p have a nontrivial 
periodicity in the direction of ra^ and are constant in the three directions perpendicular 
to it. The coordinate length of the period in ra-direction is 



Ax 



(3.5) 



with K{i) an elliptic integral of the first kind: 



K(^) 



dt 



If we employ the dimensionful coordinates 



1 l2 

r(-) 



1.31 



(3.6) 



^ /47r 

the line element related to the minimum action configuration (13. 4p reads 



(3.7) 



ds"^ = — sn^ 

c 



n 4(1 



PI 



(3.8) 



Here = (sn^(-;z)) is taken to be the average, over one period, of the squared Jacobi 
function: 



where E(z) is an elliptic integral of the second kind 

m- 



0.45 



dt\l- ^ 1.91 







1 -t2 



(3.9) 



(3.10) 



With this particular normalization of the coordinates, the average of ds"^ equals the 
standard form of the line element on flat euclidean space, 



{ds^) = S^udx^dx" ^ 



(3.11) 



even though the scale factor is rapidly oscillating on short scales. If -u = 0(1), the 
coordinate length Ax of one period is of the order of a Planck length. 

The conformal factor pertaining to Eg. (13. 8 p is depicted in Fig. ([3]). A potentially un- 
physical feature of (13. 8 p is that the metric degenerates at the zeros of the Jacobi function. 
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Figure 3. The solution (|3.4|) for x = ^, u = 12 and = S^i. 

4 The variational approach 

In this section we employ a variational technique in order to minimize the restricted action 
functional f l2.15p . Here we shall focus on the case of vanishing cosmological constant, 
A = 0, whence u = 0, and 

Slx] = lld'x{-id,xr + ^,{^y} (4.1) 

where we abbreviated 

Q ^ (72/3)^1/2 ^4 2) 

The impact of the cosmological constant will be explored later on with a different method. 



4.1 The single plane- wave ansatz 

We start by using a variational ansatz in which the modulation of the conformal factor 
has the structure of a single plane wave. Without loss of generality we may use a frame 
such that the corresponding wave vector points along the = z direction. Thus the trial 
ansatz for x depends on a single cartesian coordinate only: 

X{z) = A 1 + hcos{uz) (4.3) 

Eg. (14.31) comprises a 3 parameter family of conformal factors. They are labeled by the 
variational parameters A, h and u which we are going to adjust in such a way that (14. ip 
restricted to the trial space assumes its minimum. 
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By inserting (14. 3p into ( 14. ip . we find 

^ = Vol ■ ^ with ^ = ^1 dz[- {d,xf + ^ (^) } (4.4) 

where Vol is a (very large) four- dimensional normalization volume element, and denotes 
the action, per transverse 3- volume, averaged over one period 2ti/v in the 2;-direction. 
Explicitly, the integral for simplifies to 

dz 

SttJo 2Q'^[l + 2hcos{uz) + h?coi^{vzf) 



^ = — / dz _^^^ . _ . . . —-^ (4.5) 



- A2Q2 + ^2g2 (i, _ 2 A2Q2/^ COS {v z) + 2 A^Q^h COS (z/ zf 
-A^QV cos (z/ ^)' + A2g2/,2 cog ^ ^^2 ^2" 

Performing the integrations yields the exact result 



From this expression we obtain the following partial derivatives with respect to the two 
variational parameters h and v: 



dh 2g2(l _ ^^2)5/2 



(4.7) 



(4. 



h^A^Q'^^JY^ - K'A^Q^Vl^ + 8 K'v^ - 4 v^Vl^K' - Av"" + A u^Vl^ 

The minima of the action can now be searched by equating (14. 7p and (14. 8p to zero and 
solving for h and u. After some manipulation, the following solution is obtained: 

h = 0.5830, u = 0.4595 A Q, 5[x]/Vol = -0.008960 A^Q^. (4.9) 

This solution amounts to the non-degenerate, global minimum of ^{h, u, A) consid- 
ered a function of h and u only, with A kept fixed. Indeed, the variational principle does 
not fix the overall normalization A; as we shall see in a moment it will assume a unique 
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value once we adopt standard coordinates on the flat euclidean space which arises in the 
macroscopic limit. 

Before closing this subsection we comment on a technical issue. Due to the compli- 
cated denominator of the integrand in fl4.5p it is difficult to extend this approach to more 
complicated test functions since usually the integral cannot be performed exactly any 
longer. Let us therefore use the exactly soluble case (14. 5 p in order to test an alternative 
strategy, namely to expand the integrand in (14.50 in a power series of hnpto a given finite 
order prior to the ^-integration. This latter integration is then performed analytically on 
a polynomial expression of trigonometric functions. The expansion is found to rapidly 
converge to the values (14.91) when we increase the order of the polynomial, as it is possible 
to notice in Tab.(IT]). 





Table 1. 


Convert 


;ence 


order 


h 


i^/AQ 




0{h^) 


0.7520 


0.3908 


-0.009884 




0.6337 


0.4339 


-0.009317 




0.6016 


0.4489 


-0.009094 




0.5902 


0.4550 


-0.009010 


0{h'^) 


0.5857 


0.4579 


-0.008978 



4.2 The combination of two plane-waves 

The important question is to understand if a combination of several plane-waves can lower 
further the value of the action. In order to address this point a more general trial ground 
state needs to be considered. In particular, if we set 

x{x) = A 1 + h(^r cos{ujz) + cos(z/2; + sy)j (4-10) 

we can discuss the possibility that a combination of two plane-waves attains a lower value 
for the action by considering simultaneous variations of h, r, z/, u and s. Without loss of 
generality, the wave vector of the first plane wave has a 2;-component only, and that of 
the second lies within the z — y plane. 

In this case it is more convenient to perform the z integration after the expansion 
in powers of h of the integrand, as suggested before. The action density thus becomes 

sw/voi^^i±i^!<^;±^ (4.11) 
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where we find for the three coefficients, to order 0{h^), 

Ao = 8 z/^ + 18 h^u^ + 25 h^u^ + 90 /i + 15 /iW^ + 25 h^r'^u^ + 15 h^r^u^ 

+12 h^u'r' + 18 /i W - 8 Q'A'u^ - 8 Q'A'r^u^ + 12 /i W + 8 r V 

+120 h^rWu^ + 120 /iVwV^ + 90 /iW + 48 h^Wu^ (4.12) 

A2 = 120 /iV^wV^ + 36 h^u^ + 180 hW + 16 z/^ + 30 /iW^ + 120 h^rWu^ 
+50 /iV^ + 24 /i^i/V^ - 8 g2^2^2 ^ 4g /,2^2^2^2 ^4_^3^ 

A4 = 90 /iW^ + 8 z/^ + 18 /iV^ + 15 /iW^ + 25 h^u^ + 12 /iW^ (4.14) 

The minimization procedure involves a solution of a rather involved system of five 
nonlinear equations. It turns out that the only real solutions are the previous single plane 
wave solutions 

h = 0.7521, u = 0.3908 AQ, u = arbitrary, r = 0, s = (4.15) 
and a class of coupled plane-waves labeled by the parameter s, 



' 5QQ 135\/5 

/i = 0.4229, u=\l——- ^AQ, u; = 0.3948AQ, r = 1, s = arbitrary (4.16) 

1271(1 + s^) 

However, we find that the value of the action is lower (more negative) for the single plane 
wave solution fl4.15p than for for multi plane-wave solution f l4.16p . being — 0.0108A^(5^ in 
the first case and —0.00697^^(5^ the second case. We have checked that this conclusion 
is further reinforced by the inclusion of additional powers of h in the determination of the 
minimum up to 0(/i^^). 

4.3 The optimal trial metric 

Summarizing the above results for the case A = we can say that the variational calcu- 
lations indicates that the absolute minimum is a single nonlinear plane wave which can 
be approximated by the harmonic ansatz 

X{x) = A[l + h cos {^QA n^x'' + a)] (4.17) 

with an arbitrary unit vector and phase The constants h ^ 0.58 and = vj AQ ^ 
0.45 are universal numbers independent of /3. The overall constant A is not determined by 

^Clearly this is reminiscent of the toy model in [1] where it has been shown rigorously that the global 
minimum is assumed for a single (harmonic, in this case) plane wave. 
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the minimization condition but rather gets fixed when we adopt a normahzation condition 
for the dimensionful coordinates. Here it is natural to introduce 



x''=^/'^(2 + h^)Aipix'' 



(4.18) 



When expressed in terms of these coordinates the fine element related to fl4.17p reads 



2 + 



1 + h cos 



487r(2 + /i2)/3 



mpi n^x'^ + a 



^ fj,y dx doc 



(4.19) 



Averaging over the harmonic oscillations we find a fiat spacetime with {ds^) = S^^dx^dx^ 
again. Note that the frequency of the oscillations increases for decreasing /3. It approaches 
infinity in the limit of a pure Einstein- Hilbert action, /3 — )■ 0. Note also that, since h <1, 
the conformal factor of fl4.19p has no zeros and the metric is everywhere non-degenerate. 

5 Impact of the cosmo logical constant: 
numerical solutions 



To supplement the analysis of the previous section where we set A=0 we shall now allow 
for a nonzero cosmological constant and explore its impact on the stationary points. We 
restrict our attention to solutions depending on one cartesian coordinate only, = x, 
say. The equations f l2.2ip and f l2.25p for x and cr, respectively, are then most conveniently 
written in the form 



72/3 



.11 I ^,,3 



6 



XV"-3X(X")'-4XX'X'" + 6X"(X')' 







and likewise 



a +(a) 



7/ 

6 



a — ba [a ) 



(5.1) 



(5.2) 



Here the prime denotes a derivative with respect to x. 



5.1 Linear approximation 

Later on we shall use a numerical technique in order to find periodic solutions to these 
ordinary differential equations. Before embarking on that it is useful to analyze their 
linearization which describes small deviations from flat space, i.e. from x = 1 or cr = 0, 
respectively. 
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Figure 4. Analytical solution of eq. (j5.3p for Q = 2 and n = 5 - 10^^ (solid line). We observe 
a slow periodic "Hubble" evolution with frequency a;_, superimposed with rapid "Planckian" 
oscillations with frequency a;_|_. For comparison also the solution for Q = 10^^ and n = is 
shown (dashed line) . It represents a magnification of the single-scale "Planckian" oscillations 
for vanishing cosmological constant; it could be regarded a model for Minkowski space at 
small distances. 



Starting from Eg. (15. 2 p we expand e*^ = 1 + a + . . . and retain the derivative and 
non-derivative O(cr^) terms, as well as the 0(cr°) term. This leads to the linear equation 

a" + ^(l + 4a) + 72/3(T"" = (5.3) 
o 

Its most general solution reads 

a{x) = ~^ + ^+ cos(x c<;+) + A_ cos(x c<;_) + 5+ sin(x 0;+) + B_ sin(x u;_) (5.4) 
where A±, B± are integration constants, and 

a;± = -^\/g(3g± V9g2_24M). (5.5) 

Recall also that u = ^GK and = 1/(72/3) > 0. 

The constants uj± are real, and so cr(a;) is periodic if the cosmological constant is 
small enough, namely when > 8u/3. Otherwise they have an imaginary part which 
leads to an exponential behavior of the solution. 

For an exactly vanishing cosmological constant {u = 0) we have u+ = Q and U- = 0, 
whence a{x) is periodic with a single period oc 1//3 determined by the term. 
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When we switch on a small cosmological constant the "large" frequency cu^ = Q is not 
affected much, but cr(x) develops an additional periodicity with the "small" frequency 
U-] to leading order in m, we have U- ~ A/2n/3 oc -\/A. 

Thus in general the solution fl5.4p displays two different scales on which it varies. 
If ^ ti, — )■ and the only relevant frequency is the period of the by now familiar 
i?^-term. For large Q- values the corresponding frequency is proportional to 1/a/5, and 
since Qx oc x/{^J pi-pi) we see that dimensionful coordinate period Ax is of the order 
\/^£pi, consistent with the findings above. 

The second scale is set by the cosmological constant which can turn the oscillating 
solution (15. 4p into an exponentially one for A large enough, Q"^ < 8u/3 . 

The overall behavior of is governed by the interplay of those two relevant length 
scales. Their interpretation is particularly clear when /3 = 0(1) and A ^ mpj since then 
u/Q'^ oc GAP <^ 1. In this case the small period is of the order of the Planck length 
ip\ = VG while the large period is the "Hubble" scale oc l/y/X. Hence eq. (l5.4p describes 
a slow cosmological evolution caused by A superimposed with rapid oscillations on the 
Planck scale. An example of this situation is depicted in Fig.(jl]). This picture is confirmed 
by the numerical investigation discussed below. 

If < 8u/3 the general solution of (15. 3p is a linear combination involving growing 
and decaying exponentials which drive the solution outside of the linear regime. To see 
what happens then, we "switch off' the Planck scale oscillations and set /3 = in (15. ip of 
X = e'^ . This results in 

X" + ^X' = (5.6) 

which happens to be the one dimensional restriction of the Yamabe equation (13. 2 p dis- 
cussed earlier. (From the conceptional point of view this mathematical equivalence is 
to some extent coincidential, however. Eq. (l3.2p had the interpretation of a Bogomolny- 
like equation for the minimum action configuration for the special parameter value 7 = 
24-U/3 — 1 = which requires /3 7^ 0. Only then the solution of the Bogomolny equation 
has the lowest possible action, S = 0. In the case at hand, instead, we set /3 = so that 
(15. 6p relates to the simple Einstein-Hilbert action only.) 

Nevertheless, we know that (15.61) has the periodic solution 

Xix) = xsn {x^/u/12 X + a;i) (5.7) 

with u oc GA = £pjA. Thus, introducing x oc i-pix the first argument of the sn function 
is essentially \/ Ax. This nonlinear oscillation, with a period Ax oc 1/ \/A represents the 
generalization of the above "slow", or "cosmological" variation in the nonlinear regime. 
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Figure 5. Upper panel: numerical solution of eq. (j5.2p for u = 10^^ and Q = 2 (left) and 
Q = 0.8 (right). Lower panel: numerical solution of eq. (|5.2p for Q = 3 and u = 10^^ (left) 
and u = 10^^ (right). 

This discussion suggests the following general structure of the stationary points in 
general (and also of the absolute minimum, hopefully). If for simplicity A ^ w^pi, the 
conformal factor has a double periodicity; on small length scales it undergoes oscillations 
(with period Ax oc which are due to the higher derivative term in the action. These 

oscillations are superimposed on another type of oscillations of a much larger period set by 
the cosmological constant: Ax oc l/y/X. This feature, too, is confirmed by the numerical 
investigations to which we turn next. 
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Figure 6. Minimum action configuration of the lattice regulated action (j5.9p obtained for 
/3 = 27 and u = 0. 



5.2 Numerical solutions 

We are in particular interested in numerically determining the periodic solutions of Eq.f lS^ 
In order to achieve this goal it is convenient to follow the method presented in [ID] and 
to consider fl5.2p as the boundary value problem defined by the conditions 



being xi and X2 the initial and final limits of the integration interval, respectively. A 
problem of this type can be conveniently solved by means of the shooting method embed- 
ded in a globally convergent Newton-Rawson algorithm. In other words, the conditions at 
the outer integration extremum X2 are mapped in a functional dependence of the initial 
conditions for a' = a'{x = xi) and a'" = a"'{x = xi) in order to satisfy (15. 8p aX x = X2- 
Once convergence is achieved, it is possible to adiabatically explore the parameter space 
spanned by the two variables Q and u. 

The results are summarized in Fig. ([5]). In particular in the two upper panels the 
effect of decreasing Q (which means increasing /3) while keeping u fixed is shown, so that 
the frequency of the small scale i?^-induced fluctuations increases as Q decreases (right 
upper panel). On the other hand, if we instead keep Q fixed and change u (lower panels), 
the large scale non-linear periodicity emerges as u is decreased. 



cr(xi) = a"{xi) = 



(5.8a) 



a{x2) = (t"{x2) = 



(5.8b) 
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Figure 7. Minimum action configuration of the lattice regulated action 
/3 = 55 and u = 0. 

5.3 Lattice regulated model 



obtained for 



In more than one dimension the task of determining the general solution of fl2.2ip or 
( I2.25P is extremely hard due to their elliptic non-linear structure. A possible strategy is 
to directly minimize the lattice regulated version of the functional f l2.24p as a multivariate 
function of the field at each lattice site: 

X 

+ [e^"^^'n^(a; + e^) + a{x - e^) - 2(t{x) + (cr(x + e^) - (T{x)f) 

+ ^ 36/3((T(a; + e^) + a{x - e^) - 2a{x)){a{x + e^) + a{x - e^) - 2a{x)) 

+ ((j(a; + e^) + cr(a; - e^) - 2a{x)){(j{x + e^) - cr(x)) 
+ ((T(a; + ty) + o{x - Cy) - 2a{x)){a{x + e^) - a{x)) 
+ [a[x + e^) - a{x)){a{x + e^) - (j{x)f 



(5.9) 

In (15.91) the variable a is a dimensionless field so that the lattice cutoff is a = 1, and 
CfJ^ = 5^. For computational reasons we work in ci = 2 dimensions and we consider a 
two-dimensional lattice. 

The size of the lattices used ranges from 60 x 60 mesh points up to 120 x 120 mesh 
points, in order to check the numerical stability of the results. For actual calculations 
the Fletcher- Reeves conjugate gradient algorithm, which implements a succession of line 
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Figure 8. Minimum of the lattice regulated action (j5.9p obtained for /3 = 41 and u = 0.055. 

minimizations, turned out to be particularly convenient. After an initial search generated 
by a uniform random field distribution, the direction is chosen using the gradient of the 
action. The line minimization is thus carried out iteratively in that direction [13j and the 
algorithm works quite well when the period of the term is the dominant one. 

In Fig. ([6]) and Fig. ([7]) the field configuration corresponding to the global minimum 
of the action is displayed for m = 0, /3 = 27 and /3 = 55 respectively, corresponding to 
A = . It is reassuring to notice that it is an essentially one- dimensional object, having 
a nontrivial dependence on one coordinate only. This proves, at least for = 2, that 
the global minimum is a single non-linear plane wave, in complete agreement with our 
previous variational calculation. 

Fig.® is for A 7^ and shows instead the result of the minimization obtained for 
/3 = 41 and u = 0.005. Here, too, the global minimum action configuration consists 
of a single non-linear plane wave with a double-periodicity. Further exploration of the 
parameter space shows that generically the possible "zoo" of stationary configurations 
has this very special feature. 

6 Conclusions 

In this paper we explored the possibility that the conformal factor instability of the 
Einstein-Hilbert action is cured by an additional higher derivative invariant. As a concrete 
realization of this idea we considered A + R + gravity in 4 dimensions, for technical 
reasons restricted to the purely conformal sector. Using various techniques we found that 
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the field configurations which correspond to the global minimum of the relevant action 
functional is a family of non-linear plane waves which are labeled by a unit vector and 
a phase angle a. While the field is constant on the hyperplanes perpendicular to n^, it has 
a nontrivial modulation along the direction of n^. Approximately, when the cosmological 
constant is small, this modulation can be thought of as the superposition of two harmonic 
waves with two frequencies determined by the cosmological constant and the coefficient of 
i?^, respectively. If their scales are well separated, we are led to the picture of Planckian 
ripples due to the term which are superimposed on a much smoother solution to the 
simple K + R theory. 

This result will be a key ingredient in the analysis of the ground state of the quantized 
theory. Following the approach of [1] the logical next step will consist in a computation 
of the pertinent effective action F by means of a saddle point expansion about the global 
minimum which we found. Since this minimum is degenerate with respect to and a, 
this computation will involve an integration over the moduli space x already at 
the leading order of the semiclassical expansion, thus restoring full 0(4) invariance at the 
level of F. Since the saddle point is known only numerically, this is a rather difficult step; 
we shall came back to it elsewhere. 

The effective action functional F will in particular contain information about the 
theory's ground state and small fluctuations around it. If the semiclassical expansion is 
valid we may expect this ground state to be symmetry breaking, i.e. modulated, picking 
a specific point in the vacuum manifold 5*^ x S^. 

There is an important analogy between the underlying stabilization mechanism of 
the conformal factor instability and known examples of modulated phases in Solid State 
Physics. In particular, in either case spatially nonconstant modes condense as the result 
of the competing (negative) nearest-neighbor kinetic term and the (positive) next-to- 
nearest-neighbor term which, in gravity, stems from the four derivatives in R^. 

It would be nice to discuss possible experimental or observational signatures of this 
phenomenon, along the lines suggested in |44j, and we hope to address this issue in a 
following work. 
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